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Abstract 

We propose to search for the neutrino radiative decay by fitting a photon energy spec- 
trum of the cosmic infrared background to a sum of the photon energy spectrum from 
the neutrino radiative decay and a continuum. By comparing the present cosmic infrared 
background energy spectrum observed by AKARI and Spitzer to the photon energy spec- 
trum expected from neutrino radiative decay with a maximum hkehhood method, we 
obatined a hfetime lower hmit of 3.1 x 10^^ to 3.8 x 10^^ years at 95% confidence level for 
the third generation neutrino in the mass range between 50 meV/c^ and 150 meV/c^ 
under the present constraints by the neutrino oscillation measurements. In the left-right 
symmetric model, the minimum lifetime of is predicted to be 1.5 x 10^^ years for of 
50 meV/c^. We studied the feasibility of the observation of the neutrino radiative decay 
with a lifetime of 1.5 x 10^'' years, by measuring a continuous energy spectrum of the 
cosmic infrared background. 
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1 Introduction 



The difference between tlie mass-squares of different-generation neutrinos lias been mea- 
sured by atmosplieric neutrino oscillation experiments with Super-Kamiokande [H [2], 
neutrino beam experiments with K2K [3] and MINOS [H [5], solar neutrino oscillation 
experiments with Super-Kamiokande |6] and SNO [3 [8] and a nuclear reactor oscillation 
experiment with KamLAND [9l [10] , but the neutrino mass itself has not been measured 
yet. Detection of neutrino radiative decay enables us to measure a quantity independent 
of the difference between the mass-squares of different-generation neutrinos. Thus we 
can determine the neutrino mass itself from these two independent measurements, the 
neutrino oscillation and the neutrino radiative decay. 

As the neutrino lifetime is so long as to be much larger than the age of the universe, 
the most promising method is to observe the decay of the cosmic background neutrino. 
Detection of this radiative decay means a discovery of the cosmic background neutrino 
predicted by standard cosmology. 

The cosmic background neutrino has a temperature of 1.9K and a particle density p 
of 110 cm~^ per generation. 

Atmospheric neutrino data with Super-Kamiokande and neutrino beam data with K2K 
and MINOS gives the mass-square difference between z/3 and z/2, Am32 of (2.43 ± 0.13) x 
10~^el^^ [II] and the mixing angle of sin^26'23 > 0.92 at 90% confidence level [IT] . 

Solar neutrino data with Super-Kamiokande, SNO and the reactor neutrino data with 
KamLAND give the mass-square difference between V2 and vi of Am^2 — (7.59 + 0.19/ — 
0.21) X 10-^e1/2 [XI] and the mixing angle of sin22^i2 = 0.87 ± 0.03 [A]. 

We search for the neutrino radiative decay under these constraints on the mass-square 
differences and the mixing angles between the neutrino generations. For this search, we 
make a fit of the photon energy spectrum of Cosmic Infrared Background (GIB) to a sum 
of a neutrino decay photon energy spectrum and a continuum spectrum. 

In the left-right symmetric model, we make a feasibility study of the neutrino decay 
search in the GIB energy spectrum measurement. 



2 Neutrino Radiative Decay 

In this paper we report the analysis results on the normal hierarchy case where the z/3 is 
the heaviest neutrino and decay into z/2 or z/i as z/s — )■ z/2 ( or z/i ) + 7. 

In the z/3 rest frame, the decay photon energy Eq is related to the neutrino masses by 



Since the term m| — m\ was measured by the neutrino oscillation experiments, we 
can determine by measuring this Eq. The plot of versus m2 is shown in Fig. dJ 
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The decay photon energy Eq is given as a function of with the neutrino oscillation 

measurement constraint as shown in Fig. |2l When we consider the 7723 range between 50 

meV/c^ and 150 meV/c^ , Eq ranges from 25 meV to 8 meV in the far-infrared region. 

Taking into account the redshift of the decay photon energy, the observed photon 

energy E^ is given by 

E = ^ 
^1 1+2 ' 

where z is a redshift. As the lifetime become longer because the neutrino is moving 
away, the decay rate of neutrino R is given by R = ^^^^^^ , where r is a neutrino lifetime 
at the rest frame. Thus the decay photon flux per unit solid angle is given by 

dSdtdn 47rT(l+2)"'' ' 

where r is a distance between a neutrino decay point and a telescope with a detection 
area dS, and p is a density of the cosmic background neutrino of 110 cm~'^. Assuming a 
flat universe, dr is related with dz by 

dr = -^^[(l + zfilM + ^Aj-^^-^dz, 

where Hq is a Hubble constant of ( 74 ± 4 ) km s~^ Mpc~^, and Qm and are the 
matter density and the cosmological constant which were measured to be [11] 
= 0.26 ± 0.02, 
= 0.74 ±0.03, 

riM + nA = 1.006 ± 0.006. 

Thus under the assumption of a flat universe with Qm + = 1, we obtain the 

dN^ 
dSdt 



following distribution of the decay photon flux per unit solid angle and unit energy : 



A 



dSdtdndE^ 4tttHoE^ E. 

If we use the energy flux per unit solid angle / deflned by 

T _ d{NryXEj) 

— dSdtdfl ' 

we obtain the energy flux per unit solid angle and unit energy given by 

dE-, - A-ktHo \\e~,I ^ ^A/ + ^ i Aj 

Using E^ = hu , 

This spectrum is smeared by the neutrino motion at 1.9K, but this effect is negligibly 
small. The photon energy spreads due to this neutrino motion are 0.6% and 1.5% in rms 
for the photon energies of 25 meV and 10 meV, respectively. 



3 Neutrino Lifetime Limit 

The Cosmic Infrared Background (CIB) continuum in photon energy spectrum is the 
most serious foreground against this neutrino radiative decay signal because the sharp 
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edge of the signal spectrum is located between 8 meV and 25 meV in our search region. 
There have been studies where the neutrino lifetime limit was obtained using the CIB 
data [121 IE] • However they did not use the photon energy spectrum expected from the 
neutrino radiative decay to estimate the neutrino lifetime limit but used only the total 
decay rate of neutrino. They did not use the AKARI CIB data which were not available 
at that time. 

We perform a statistical test of the CIB data measured by COBE [T^ [T5] and AKARI 
[16] to estimate the neutrino lifetime limit. Since we do not have a sharp edge at high 
energy end of the CIB energy spectrum, we set a lower limit of the heaviest neutrino 
lifetime. The CIB data measured by COBE and AKARI were directly compared with the 
photon energy spectrum expected from the neutrino radiative decay neglecting the CIB 
continuum with a maximum likelihood method. In this comparison, we neglected the CIB 
continuum in order to estimate the neutrino lifetime limit in the worst case which gives 
the lowest lifetime limit. The curve with a maximum likelihood is shown together with 
the likelihood as a function of lifetime for = 50 meV/c^ and m2 = 10 meV/c^ in 
FigO The arrow point in this figure gives us a lower limit of lifetime at 95% confidence 
level. We performed this procedure for various masses to obtain a lower limit of the 
1^3 lifetime at 95% confidence level as a function of mass as shown in Fig. HI They 
range from 1.7 x 10^^ to 2.4 x 10^^ years in the mass range between 50 meV/c^ and 
150 meV/c^. 

Recent deep galaxy surveys with infrared satellites, AKARI [16], Spitzer [17] and 
Hershel [18] have revealed that more than a half of the CIB energy originates from external 
galaxies. So we performed the analysis to obtain stronger constraint on the neutrino 
lifetime taking the integrated flux of galaxies into account. We subtracted the contribution 
of the distant galaxies as source points [17] from the CIB measured by AKARI as shown 
in Table [H Since there were no measured values of the contributions at wavelengths of 
60, 90 and 140 /xm, we interpolated the measured values at wavelengths of 24, 70 and 
160 fim. Thus we obatined the CIB after subtracting the contribution of distant galaxies. 
These corrected CIB data give a constraint on the neutrino lifetime as the most stringent 
lower limits. 

With the AKARI CIB data after subtracting the contributions of distant galaxies 
as point sources, we performed a statistical test of the spectrum to obtain the neutrino 
lifetime limit, motivated by the fact that such corrected AKARI CIB data has much less 
contributions of distant galaxies than the COBE CIB data. The corrected AKARI CIB 
data were compared to the photon energy spectrum expected from the neutrino radiative 
decay with a maximum likelihood method. The curve with a maximum likelihood is 
shown together with the likelihood as a function of z/3 lifetime for rris = 50 meV/c^ and 
m2 = 10 meV/c^ in Fig [51 An arrow point in this figure gives us a lower limit of lifetime 
at 95% confidence level. We performed this procedure for various u-^ masses to obtain a 
lower limit of the lifetime at 95% confidence level as a function of z/3 mass as shown in 
Fig. O They range from 3.1 xlO^^ to 3.8 x 10^^ years in the z/3 mass range between 50 
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meV/c^ and 150 meV/c^. 



4 Neutrino Lifetime in the Left-Right Symmetric SU(2)/^x 
SU(2)i?x U(l) Model 

In the standard model, the heaviest neutrino hfetime is predicted to be 10^^ year for z/3 
with a mass of 50 meV/c^[T9l |20l |2T1 122]. It is too long to be measured by the present 
method. 

In the left-right symmetric model, the lifetime is predicted to be much shorter than 
in the standard model. The left-right symmetric SU(2)iX SU(2)j:jX U(l) model has two 
charged weak bosons, or the left-handed weak boson Wl and the right-handed weak boson 
Wfi which are mixed with a mixing angle ( into two mass eigenstates Wi and W2 as follows 

Wl = Wlcos( - WfisinC 
W2 = WLsm( + VTrcosC 

In this model, the lifetime of 1^3 is given by 

= ^i'^flUs.msnU^l + ml)^^{l + + 4m?(l - sin^2C 

where a is a fine structure constant, is a Fermi coupling constant, m,-, Mwi and 
Mw2 are masses of r, Wi and W2, respectively [211 122]. Uij is the (i, j)-th element of the 
Maki-Nakagawa-Sakata mixing matrix [23] and we took |?732| = 1/a/2 and {U^sl = l/\/2. 
The present mass limit of the right-handed weak boson Wr is Mr > 0.715TeV assuming 
C = 0. The upper limit of the mixing angle between Wl and Wr is sm( < 0.013 [H]. The 
lifetime is shown as a function of rris with of 0.715TeV, sin( of 0.013 and Amgg of 
(2.43 ± 0.13) X IQ-^eV^ in Fig. [71 

For 777,3 of 50 meV/c^, the lifetime is calculated to be 1.5 x 10^^ years under the above 
conditions in left-right symmetric model. The present lifetime limit is shorter than this 
prediction by a factor of 2 xlO~^. 

We can improve the sensitivity of the search for neutrino radiative decay by measuring 
the continuous energy spectrum of the Cosmic Infrared Background (CIB). We fit the 
continuous spectrum with a sum of the CIB continuum and the neutrino decay photon 
spectrum to search for a sharp edge of the neutrino decay photon spectrum. 

To estimate the energy resolution requirement, we performed the simulation study of 
the neutrino radiative decay assuming the following conditions: 

• The present CIB continuum is dominated by other sources than the neutrino ra- 
diative decay. This CIB spectrum is represented by a Planck distribution plus a 
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quadratic function. 

• m3 = 50 meV/c^ and m2 = 10 meV/c^. 

• rs = 1.5 X 10^^ years. 

We simulated the CIB detection experiments with a 20cm-diameter telescope with a 
viewing angle of 0.1 degrees and energy resolutions of to 5 % by 1 % step. With 100% 
detection efficiency and 10- hour data taking, we have the photon energy spectrum for the 
CIB and the neutrino radiative decay as shown in Fig. [HI The CIB continuum was fitted 
to a function of the Planck distribution plus a quadratic function. Then we calculated the 
energy derivative of the photon energy spectrum for a sum of the CIB continuum and the 
neutrino radiative decay as shown in Fig. EJ In this negative energy derivative plot, the 
sharp edge of the energy spectrum is identified as a clean peak. The excess is 6.7a and 
we can see the clear signal peak in the distribution if the energy resolution is less than 
2%. 

We search for the neutrino radiative decay in the range between 50 meV/c^ and 
150 meV/c^ which corresponds to the photon energy range between 25 meV and 8 meV 
with this experiment. Thus the photon energy range to be measured is between 5 meV 
( A = 250 fim ) and 35 meV ( A = 35 /im). In this photon energy range, the frequency 
of photon coming in this telescope is expected to be around 5 MHz. As we will use 400 
pixels as a photon detector, the photon rate per pixel is around 12 kHz/pixel. 

By this simulation study, we found that the required energy resolution is less than 
2% at 25 meV. With this resolution, we estimated 5a observation lifetime with 10-hour 
running of this telescope from the differential CIB energy spectrum including the radiative 
decay of the cosmic background neutrino for the range between 50 meV/c^ and 150 
meV/c^ with Amgg of 2.43 x lO'^e^^. The 5a observation lifetime ranges from 2.0 x 10^'' 
to 3.0 X 10^^ years as shown in Fig. [71 Spectral emission features of the dusts in galaxies 
at redshifts of 1, when a large fraction of the CIB energy was generated ^7\, may produce 
a structure in the spectral energy distribution (SED) of CIB, similar to the neutrino 
decay photon SED. However, SEDs of such distant galaxies in the energy range where the 
neutrino decay photon is expected to be observed, have not been well explored. Further 
study of the SED of distant galaxies with future, large-aperture, infrared telescope such 
as JAXA's SPICA (Space Infrared Telescope for Cosmology and Astrophysics) mission 
[23] is important to mitigate the spectral contamination of the neutrino decay photon. 

The present measured CIB spectrum includes the point-sources of distant galaxies and 
the zodiacal light foreground ambiguity. By using small pixels with small viewing angle 
such as 0.005 degrees, we will be able to distinguish the distant galaxy point sources from 
the CIB spectrum. In the future project EXZIT ( Exo-Zodiacal Infrared Telescope ) [25] 
which will observe the CIB outside of Jupiter orbit, we will decrease the effect of the 
zodiacal emission significantly, and will have much less ambiguity of the zodiacal light 
foreground. 
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5 Conclusion 



We propose to search for the neutrino radiative decay by fitting a photon energy spec- 
trum of the cosmic infrared background to a sum of the photon energy spectrum from 
the neutrino radiative decay and a continuum. By comparing the present cosmic infrared 
background energy spectrum observed by AKARI and Spitzer to the photon energy spec- 
trum expected from neutrino radiative decay with a maximum hkehhood method, we 
obatined a lower hmit of netrino hfetime of 3.1 xlO^^ to 3.8 xlO^^ years at 95 % confi- 
dence level for in the mass range between 50 meV/c^ and 150 meV/c^ under the present 
constraints by neutrino oscillation measurements. 

In the left-right symmetric model, the lifetime is predicted to be much shorter than in 
the standard model. The present mass limit of the right-handed weak boson Wr is rriR > 
0.715 TeV and the mixing angle between Wl and Wn is sin^ < 0.013 . In this model, the 
minimum lifetime of is predicted to be 1.5 x lO^'' years for rris of 50 meV/c^. 

The present lifetime limit is shorter than this prediction by a factor of 2 x 10~^. 
By measuring the continuous energy spectrum of the cosmic infrared background in the 
photon energy region between 5 meV (A = 250 /xm) and 35 meV (A = 35 /im), we can 
expect 5a observation of the radiative decay of neutrino with a lifetime of 1.5 x 10^^ years. 
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m2 (eV) 

Figure 1: Relation between ms and The solid band shows the la constraint by 
the neutrino oscillation measurements (Amgg of (2.43 ± 0.13) x lO^^eV^ ). Three curves 
correspond to various neutrino decay photon energies of 25meV (dashed), 20meV (dotted) 
and 15meV(dot-dashed). (Color online) 
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Figure 2: Relation between Eq and m^. The solid band shows the la constraint by the 
neutrino oscillation measurements (Am|2 of (2.43 ± 0.13) x lO^^e^^ ). 
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Figure 3: Top plot shows the CIB data measured by COBE (dark square) and AKARI 
(dark circle) fitted to the photon energy spectrum expected from the neutrino radiative 
decay for 7723 = 50 meV/c^ and m2 = 10 meV/c^ with a maximum likelihood method. 
The curve shows the best fit. Bottom plot shows the likelihood as a function of the 1^3 
lifetime. An arrow points to a lower limit of the neutrino lifetime at 95 % confidence level. 
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Figure 4: Lower limits of the neutrino lifetime at 95 % confidence level as a function of 
rris obtained with the CIB data measured by COBE and AKARI. 
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Figure 5: Top plot shows the AKARI CIB data after subtracting the contribution 
of distant galaxies as point sources (dark circle) fitted to the photon energy spectrum 
expected from the neutrino radiative decay for = 50 meV/c^ and m2 = 10 meV/c^ 
with a maximum likelihood method. The curve shows the best fit. Bottom plot shows 
the likelihood as a function of z/3 lifetime. An arrow points to a lower limit of the neutrino 
lifetime at 95 % confidence level. 
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Figure 6: Lower limits of the neutrino lifetime at 95 % confidence level as a function of 
rrts obtained with the CIB data measured by AKARI after subtracting the contribution 
of distant galaxies. 
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Figure 7: Expected lifetime minimum in the left-right symmetric model as a function 
of 7713. The solid band shows the la constraint by the neutrino oscillation measurements 
(Am|2 of (2.43 ±0.13) x 10~^eV^ ). The sensitivity of 5 a observation with the proposed 
measurement is shown by dark circles. 
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Figure 8: The top plot shows the photon energy spectra of the CIB and the neutrino 
radiative decay with various energy resolutions from % to 5 % by 1 % step. The 
CIB continuum is fitted to a sum of a function of the Planck distribution and a quadratic 
function. The bottom plot is the negative energy derivative of the photon energy spectrum 
for a sum of the CIB continuum and the neutrino radiative decay for energy resolutions 
of 0, 1, 2 and 3 %. (Color online) 
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Wavelength 


CIB" 


contribution of 


CIB 


CIB 






distant galaxies^ 


after subtarction'^ 


after subtraction'^ 


(/xm) 


(MJy sr-i) 


(MJy sr-i) 


(MJy sr^i) 


(nW sr~^) 


24 




0.017 ± 0.003 






60 


0.27 ± 0.20 


( 0.112 ± 0.018 ) 


0.16 ± 0.20 


7.3 ± 9.2 


70 




0.138 ± 0.023 






90 


0.67 ± 0.15 


( 0.234 ± 0.042 ) 


0.44 ± 0.16 


14.4 ± 5.1 


140 


0.94 ± 0.17 


( 0.475 ± 0.095 ) 


0.47 ± 0.20 


9.8 ± 4.1 


160 


0.73 ± 0.21 


0.571 ± 0.120 


0.16 ± 0.24 


3.0 ± 5.0 



Table 1: The AKARI CIB data before and after subtracting the contribution of distant 
galaxies. 

" The CIB measured by AKARI in unit of MJy sr^^ 

^ The contribution of distant galaxies to the CIB measured by Spitzer in unit of MJy sr~^. 
The numbers in parentheses were obtained by a linear interpolation. 
^ The CIB measured by AKARI after subtracting the contribution of distant galaxies in 
unit of MJy sr^^. 

The CIB measured by AKARI after subtracting the contribution of distant galaxies in 
unit of nW m~^ sr~^. 
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